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Abstract 

Based on previously published experimental observations and mathematical mod- 
els for Hesl, p53 and NF-kB gene expression, we improve these models through a 
distributed delay formulation of the time lag between transcription factor binding 
and mRNA production. This description of natural variability for delays introduces 
a transition from a stable steady state to limit cycle oscillations and then a second 
transition back to a stable steady state which has not been observed in previously 
published models. On the basis of our results and following recent discussions about 
the role of delay-induced oscillations in gene transcription we establish the hypothe- 
sis that the period of the delay-induced cyclic changes should be characterized by an 
upper bound so that it cannot be greater then the period of fundamental biological 
cycles. We demonstrate our approach for two models. The first model describes 
Hesl autorepression with equations for Hesl mRNA production and Hesl protein 
translation. The second model describes Hesl repression by the protein complex 
Gro/TLEl/Hesl, where Gro/TLEl is activated by Hesl phosphorylation. 

Keywords: gene expression with negative feedback, distributed time delay, limit 
cycle oscillations 
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Introduction 



The three proteins Hesl, p53, and NF-kB are transcriptionally regulated by short 
negative feedback loops and actual experiments have revealed that their genes show 
an oscillatory expression 0, 0, Q • In the vertebrate segmentation clock several 
cycling genes are involved in an oscillatory mechanism driving somite segmentation 
£ ] . Mathematical models have been established for the regulation of these proteins 
I H H H 11 II E3, EH, Q where the authors of 0, EJ have introduced 
a discrete time delay to describe the time lag between transcription factor binding 
and gene transcription. Analytical and numerical results have shown that these 
models, which display damped oscillations for small positive delay times, can pass 
a Hopf bifurcation at a critical delay time where the stable steady state becomes 
unstable and a stable limit cycle oscillator emerges EJ . 

Naturally, the question arises what the role of such oscillations in gene tran- 
scription could be 0, Q, 0, 0]? Due to the fact that the oscillation period is 
proportional to the delay time it may be possible that oscillations encode informa- 
tion. Single cell experiments have shown that the dynamics of NF-kB dependent 
transcription correlates with the duration of p65 oscillations 0, ^| . In other single 
cell experiments p53 has been expressed in a series of discrete pulses in response 
to DNA damage and the average number of pulses over all cells increased with the 
damage 0. 

We introduce in our work a new model for gene transcription networks with a 
negative feedback loop where the time lag between transcription factor binding and 
mRNA production is described, in line with the natural variability, by a distributed 
time delay. With an increasing mean delay time this model can show a Hopf bi- 
furcation leading to a transition from a stable steady state to an unstable steady 
state surrounded by a stable limit cycle oscillator and then a reverse Hopf bifurca- 
tion leading to a transition from a stable limit cycle oscillator to a stable steady 
state. Such a property has not been observed in previously published mathematical 
models for gene transcription networks. Following recent discussions about infor- 
mation encoding in delay-induced oscillations our results could be interpreted that 
the period of the delay-induced cyclic changes cannot be greater then the period of 
fundamental biological oscillations, say the cell cycle or the circadian clock. There- 
fore for a large delay time a reverse Hopf bifurcation should lead back to a stable 
steady state. This is demonstrated in our model. 

The outline of the paper is as follows. The first section introduces differential 
equations with discrete and distributed time delays. In our model the distributed 
time delay is defined by a kernel of the Gamma function. The following two sections 
validate the approach with two models from the literature. The first model describes 
Hesl autorepression with equations for Hesl mRNA production and Hesl protein 
translation. The second model describes Hesl repression by the protein complex 
Gro/TLEl/Hesl, where Gro/TLEl is activated by Hesl phosphorylation. Our an- 
alytical and numerical results are discussed in relation to experimental data. This 
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leads us to an investigation into thresholds for damped and sustained oscillations. 



Discrete and distributed time delays 



An appropriate mathematical representation to describe a time lag between an 
action and a reaction are delay differential equations 17]. If a variable influences 
the present state of the system at a fixed time point in the past, t — r, where r is 
the delay time, then the time lag can be modelled as a discrete delay 



dx 
~dt 



F(x,x(t - t)) . 



A more realistic description of time lag effects in biological systems is to assume that 
the past influences the present state over an interval of time. Natural variability is 
weighted by a probability density g, leading to a distributed delay 



dx 
~dt 



F 



x(t — t)q(t)cIt 



(1) 



Analyzing a discrete delay differential equation through a Taylor series expansion 
of x(t — r) leads to an infinite dimensional system of ordinary differential equations. 
Approximations by a finite number of ordinary differential equations have been 
studied in 2> lU 2p| . In case of a distributed delay, we discuss below an example 
where a finite set of ordinary differential equations can exactly be derived. In either 
case, the increased dimension can lead to a stabilization or destabilization of the 
dynamics in feedback models. 

Delay differential equations have been extensively applied to biological and engi- 
neering processes. This includes gene transcription B , la, 1^. [lcL 21, 22] j the nucleo- 
cytoplasmic translocation of proteins in eucaryotesj2J, |2J|; the hatching period 
and the incubation time in population models 25, 2q|. In control theory the delay 



between the observation and the control of a variable has been considered |27l . l2£ 

A frequently used model for distributed time delays in biological applications is 
to choose for the probability density g in Eq. Q the kernel of the Gamma function 



(p-1)! 



(2) 



which has a maximum at r = p/q and is zero at r = and for r — > oo. The mean 
delay time is f = (p + l)/q. This leads us to the following model 



F x 



x(t - T)g p Jr)dT 



dx 
~dt 



This model has the advantage that using the linear chain trick [1£|] we can trans- 
form the delay differential equations into a finite number of ordinary differential 
equations. 
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Authors that used this idea have found that distributed delay differential equa- 
tions lead to a larger range of stability than discrete delay differential equations 
26, 29|, |30j. This is important in ecological models where the destabilization of the 
dynamics can lead to an extinction of species [26] or in physiological models where 
the disease is characterized by a stability change of a variable |29( . 



Hesl autorepression (Model I) 

A mathematical model for Hesl protein autorepression with a discrete time delay 
has been introduced in 0,0,0,0]. The discrete time delay describes the time lag 
between Hesl protein binding to the regulatory DNA and the production of Hesl 
mRNA. For a critical delay time this model passes a Hopf bifurcation where the 
stable steady state is left and the system moves to a stable limit cycle 0, El, El, fioj l . 
We here modify this model by replacing in Eqs. (3.1) and (3.2) of Ref. [lfj the 
discrete time delay by a distributed time delay with Gamma kernel of Eq. ©, 
leading to 

d b ■ k h 

—mRNA = —r T~Fnci ; ; n/ — , , — a ■ mRNA 

dt k h + (/ °° Hesl (t - t) ■ g p q (r)dr) h (3) 

1 

—Hesl = b ■ mRNA - a ■ Hesl . 
dt 

The first equation describes the temporal change of the Hesl mRNA concentra- 
tion, which depends on the autorepression by the Hesl protein. The second equation 
describes the translation of the mRNA. We have implemented the distributed delay 
such that the averaging is performed over the transcription factor concentration re- 
flecting the varying binding to their sites and the dissociation. In the linear case the 
smaller translation delay can be shifted into the transcription term without chang- 
ing the dynamics 0]. Without loss of generality, we simplify Eqs. (3.1) and (3.2) 
of Ref. [lfj by choosing the same rate constants for transcription and translation 
and also the same rate constant for degradations (cf. Table We have chosen the 
parameter for the degradation rates ten-fold higher than in 0] so that the delay 
times which induce limit cycle oscillations are in the range of the temporal duration 
of gene transcription (10-40 min). 

The first step in analyzing the properties of differential equations is to look for 
steady states and to study their stability. For k = (b/a) h ■a/(b — a) the steady state 
of Eqs. © is mRNA* = 1 and Hesl* = b/a which is stable for the case without 
delay |lfj |. Linear stability analysis of Eqs. (jHJ at steady state leads to the following 
eigenvalue equation 

(A + a) 2 + — (b - a) gP{r) ■ e'^dr = . 
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The integral can be analytically solved, resulting in 



(A + a) 2 + ^(6-a) 
b 



= 0. 



The eigenvalue equation is independent of how the averaging of the distributed 
time delay in Eqs. © is performed. The eigenvalue equation is the same regardless 
whether h is outside the integral as in Eqs. © or whether h is an exponent of 
Hesl{t — t). The eigenvalue equation is also the same if we average over the 
transcription rate, i. e. , 



But it is different to the eigenvalue equation for the related model with a discrete 
time delay, Eq. (A13) in Thus considering natural variability in a delay differen- 
tial equation model for Hes-protein autorepression can lead to different steady state 
dynamics, but the steady state dynamics is independent of the source of natural 
variability. 

The next step is to analyze how the delay time influences the stability of the 
steady state by identifying critical mean delay times r = (p + l)/q, which form 
Hopf bifurcation points. First we apply the graphical method from Ref. ^| to get 
an immediate overview of putative Hopf bifurcations. We transform the eigenvalue 
equation with A = a + iu> and a = to 



Now the delay is separated into the delay curve on the right-hand-side (r.h.s.) of 
Eq. (j3J) and the ratio curve on the left-hand-side (l.h.s.). Both curves are plotted 
as a function of to in the complex plane in Fig. ^ together with the unit circle as the 
discrete delay curve. The shape of the delay curve does not depend on the scaling 
factor q. Intersections between a delay curve and a ratio curve are putative Hopf 
bifurcations. For several parameter pairs (p, h) two intersections occur indicating 
putative steady state - limit cycle - steady state transitions. For low values of p 
the system remains stable for low h or can pass through a steady state - limit cycle 
- steady state transition for higher h. For intermediate values of p three scenarios 
are possible: stable steady state for low h, steady state - limit cycle - steady state 
transitions for intermediate h and steady state - limit cycle transition for high h. 
For p — > oo the delay curve converges to the discrete delay curve. 
For p = 1 and p = 2, Eq. Q can be analytically solved for q 




b(iuj + a) 2 



(4) 



a 2 h(b - a) (iu + q)P 





with 



— (6 -a) and c p=2 = -r (b - a) . 
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The critical mean delay times r m in and r max for stability changes can be calculated 
from q± 

p+l _ _ p+l 

Tmin — and T max — . 

q+ q- 

The parameter q in Eqs. is chosen as the bifurcation parameter related to f. 
The results for f m i n and f max are summarized in Table HI This is compared to the 
results for the discrete delay according to Eq. (3.4) of |lOj. The delay times of the 
two putative Hopf bifurcation points f m in and f max are in the range of the duration 
of gene transcription. In addition, the results for a = 0.03 (parameter value in 
0, liol p are shown in Table |2J For a = 0.03 only the delay time of the first Hopf 
bifurcation point f m in is in the range of the duration of gene transcription. 

Comparing the distributed delay with the discrete delay shows that the dis- 
tributed delay destabilizes at a higher critical mean delay time than the discrete 
delay. This is in agreement with the results in Refs. 2(], 2JJ, 3(|, where the authors 
have found that the distributed delay leads to a larger range of stability than with 
a discrete delay. 

The following last step in our analysis is to check whether a changes its sign at 
the critical values f m % n and f max . For most cases this is difficult or impossible to 
calculate analytically, which is why we used numerical simulations to check whether 
the two predicted Hopf bifurcations are passed. We transform Eqs. (j3J) with p = 2 
using the linear chain trick [18( to obtain 

d 

—x = q ■ (Hesl - x ) 

^xi = q ■ (x - x{) 

d b-k h ( 6 ) 

—mRNA = — p - a ■ mRNA 

dt k h + x\ 

—Hesl = b ■ mRNA - a ■ Hesl . 
dt 

Equations (jHJ) represent a monotone cyclic feedback system. Several authors 
have focused on the existence of periodic solutions in such kind of equations [3lL 
E3- I33I 13]. According to the Poincare-Bendixson theorem for monotone cyclic 
feedback systems a stable solution of Eqs. © is either a stable steady state or a 
stable limit cycle oscillator Results from numerical simulations of Eqs. (|H]l 

with h = 6 are presented in Fig. [5J The time course of the mRNA and Hesl protein 
concentrations for delay times below and above the critical mean delay times show 
that the bifurcation at Tmin leads to a transition from a stable steady state to an 
unstable steady state. The latter is surrounded by a stable limit cycle oscillator. 
The reverse bifurcation at f max leads back to a stable steady state. Table 01 shows 
the frequency and the amplitude of the limit cycle oscillations as a function of the 
mean delay time. The frequency is inverse proportional to the delay time while the 
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amplitude has a maximum between the bifurcation points. Information can only be 
encoded in the frequency of the oscillations. 

It follows that our model of Hesl autorepression with a distributed delay has a 
maximal delay time for delay-induced limit cycle oscillations. This result is different 
to the discrete delay models [ESSE] which either show a steady state - limit cycle 
transition or no bifurcation. This can straightforwardly be shown for a generalized 
eigenvalue equation for these models, which reads (A — A) • (B — A) — C • e~ Ar = 0. 
The solution for purely imaginary eigenvalues A = iui is w 2 = — {A 2 + B 2 )/2 ± 
y/{{A 2 + B 2 ) 2 /4 - A 2 B 2 + C 2 ) and r = l/u ■ aiccos{{AB - uj 2 )/C). Multiple solu- 
tions for t due to the periodicity of the angular function do not change stability 
Thus these two-dimensional negative feedback models for gene transcription with 
discrete time delay either show no Hopf bifurcation or they can show only one Hopf 
bifurcation but they can never show a Hopf bifurcation and a reverse Hopf bifur- 
cation with increasing delay time. In contrast, a damped harmonic oscillator as an 
example from physics modelled with a discrete time delay can also show alternating 
stability and instability with increasing delay time 



The steady state properties of our model are independent on the way the aver- 
aging is performed in the distributed time delay. But they can depend on the kernel 
of the delay distribution. An alternative approach for a distributed time delay has 
been studied in Q • The averaging is performed over the transcription rate but with 
a uniform kernel of the distribution on a finite support. In contrast to our results 
the results in Q are indistinguishable from the discrete delay case. From our point 
of view, the Gamma kernel is a more realistic description for the probability den- 
sity of a distributed time delay, because it has a maximum and approaches zero for 
r — > and r — > oo. 

The infinite support of the delay distribution in our model does not lead to an 
interference of the delay induced oscillations with fundamental biological oscillations 
for a broad parameter range because with increasing delay time the kernel of the 
Gamma function exponentially converges to zero. We calculate a boundary max 
such that the contribution of the tail of the delay distribution can be neglected: 



We consider only the expression with the highest contribution of the tail in this 
sum. Assuming Hesl = 10 as an upper boundary for the protein concentration. 
Lower boundaries for the parameter max for which the inequality 



holds are max = 176min for h = 6, r = 16min and max = 672min for h = 
10, f = 48min (critical delay times of reverse Hopf bifurcations). They are smaller 
for example then a typical time interval between two cell divisions in mammals of 
24 hours or the circadian clock period. 





h ■ 
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Gro/TLEl mediated repression of Hesl (Model II) 



To see whether other models with a negative feedback loop can show a steady state 
- limit cycle - steady state transition, we consider what is referred to as Model C 
in ^jj. Here Hesl repression is realized by the protein complex Gro/TLEl/Hesl, 
denoted GroH, where Gro/TLEl is activated by Hesl phosphorylation. We replace 
in Eqs. (4.1) and (4.2) of Ref. [10] the discrete delay by a distributed delay with the 
Gamma kernel 

± R 1= b_kH H 1 

dt nesl kh + {^<f q {r).GroH{t-T)dr)h 

±GroH= 1 b h H : slh ih -a.GroH. 
dt l h + Hesl h 

For k h = a/(b — a) and l h = (b — a)/a the steady state is Hesl* = 1 and GroH* = 1. 
Linear stability analysis leads to the following eigenvalue equation at the steady state 

, 9 a 2 h 2 ,, , 9 q p 

° = < A+a > +—<'-"> Vt^ 

Again we apply the graphical method described above to get an overview of putative 
Hopf bifurcations or reverse Hopf bifurcations. Transformation of the eigenvalue 
equation with A = a + iu and a = leads to 

b 2 (iu + a) 2 q p 

a 2 h 2 (b-a) 2 (iu; + q)P {> 

The ratio curve (l.h.s.) and the delay curve (r.h.s.) of Eq. (jBJ) are shown in Figured! 
Several parameter pairs (p, h) show two putative Hopf bifurcations indicating stable 
steady state - limit cycle - stable steady state transitions. 

The critical parameter q can be analytically calculated from Eq. (JSJ) for p = 1 
and p = 2 



with 



q ± = a[c p -l±J(c p -lf-l) (9) 



h 2 h 2 
c p= i = -^{b - a) 2 and c p=2 = 7^(6 - o^ 2 



Ab 2 " 1 p " z 2b 2 

The results for the critical mean delay times r mi - n and f max are summarized in Table 
together with the result for the discrete delay (according to Eq. (4.4) in ^3]). The 
delay times of the two putative Hopf bifurcation points r m j n and f max are for some 
(p, h) pairs close to the range of the duration of gene transcription. The results for 
a = 0.03 are shown in Tabled For a = 0.03 only the delay time of the first Hopf 
bifurcation point f m in is close to the range of the duration for gene transcription. 

Comparing the results for the distributed delay with those for the discrete delay 
shows that the distributed delay destabilizes at a larger critical mean delay time 
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than the discrete delay. This is in agreement with the results for Model I. Comparing 
now the critical delay times of Model I and II shows that Model II destabilizes at 
a lower critical delay time f m i n and stabilizes at a higher critical delay time f max 
than Model I. This is related to the stronger curvature of the ratio curve of Model 
II, which comes from cooperative repression together with cooperative activation 
whereas in Model I only cooperative repression is included. 

For numerical simulations, we use the linear chain trick to derive for p = 2 the 
ordinary differential equations 

^x = q ■ (GroH - x ) 

^xi = q ■ (x - xi) 

d b-k h ( 10 ) 

—Hesl = —, r — a ■ Hesl 

dt k h + x\ 

d b-Hesl h 

—GroH = —r r — a ■ GroH , 

dt l h + Hesl h 

The simulations confirm the two Hopf bifurcations as shown in Figure HI The 
bifurcation at f m i n leads to a transition from a stable steady state to an unstable 
steady state, surrounded by a stable limit cycle oscillator. The reverse bifurcation 
at T'max leads to a transition back to a stable steady state. 



Comparison with experimental results 

In signaling and damage repair a fast response is desirable, which suggests that 
the initial reaction period is most important. This means that both, sustained and 
damped oscillations, are relevant for information encoding. The authors of Ref. 0] 
have measured discrete pulses for protein concentrations with fixed amplitude and 
duration, which do not depend on the amount of the stimulus. The experimental 
results in Ref. [f| show only few periods of the oscillations, making it difficult to 
decide whether one deals with sustained or damped oscillations. The same may be 
said for the results in '2] . A further complication arises from the difference between 
single-cell measurements, compared to averaging measurements. From a modelling 
perspective, sustained oscillations that arise after passing a Hopf bifurcation point, 
have a clear delay time threshold for frequency encoding signal transduction. 

We now discuss possible thresholds for frequency encoding in damped oscil- 
lations. The mathematical models of Eqs. (J3J) and (J7J) show damped oscillations 
for large ranges of parameters without passing a Hopf bifurcation at a critical de- 
lay time. For other parameter ranges they approach a constant steady state. If 
damped oscillations are indeed relevant then the question after a threshold arises. 
A threshold is important to filter out noise which is an inherent property of biologi- 
cal signaling pathways. For damped oscillations a threshold could be the amplitude, 
which is however more sensitive to noise than the frequency of oscillations. 
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An alternative threshold for damped oscillation could be the critical delay time at 
which for the first time at least two eigenvalues switch from real numbers to complex 
numbers; while the real parts of all eigenvalues remain negative. Such a threshold 
can be calculated analytically for a one-dimensional Hesl protein autorepression 
model (Model III) 

d b-k h 

—Hesl = — : r-z—. — jT— ; — n a, • Hesl (11) 

dt k h + (/ °° gfftr) ■ Hesl (t - T)dr) h y ' 

Linear stability analysis of Eq. Qll|) with p = 1 leads to the following eigenvalue 
equation 



resulting in 



= A + a + c 



^1/2 




2 v V 4 

From the zeros of the term under the square root 



9i /2 = a + 2 • c ± ^({a + 2 • c) 2 - a 2 ) (12) 

we calculate the mean delay times f = (p + l)/q for transitions between real and 
complex eigenvalues. For f below and above the transition values numerical sim- 
ulations of exactly derived ordinary differential equations from Eq. (|11[) show no 
qualitative difference in the dynamics of the Hesl protein concentrations, as shown 
in Figure On the other hand, numerical simulations of Model I with parameters 
h = 6, p = 2 and a = 0.4, display for intermediate delay times damped oscillations 
with higher natural frequency, while for no and large delay times the oscillations 
have a lower natural frequency, as shown in Figure H3 The eigenvalues related to 
these parameters are complex for all delay times r > but no Hopf bifurcation is 
passed. The results shown in Figures |S] and lead us to the conclusion that that 
the first switch of two eigenvalues from real numbers to complex numbers is not an 
appropriate thereshold. 

Based on the results shown in Figures El and El we suggest as a quantitative 
measure for information encoding in the frequency of damped oscillations the ratio 
of frequency and damping defined by 

< = z. 

7 

The values for u and 7 are obtained from curve fits for the data in Figures 03 and 
using the following equation 

y = a ■ e~ 7 * • cos(a;t + tp) . 



This is motivated by the concept of a damping ratio in control engineering 27]. The 
results for £ are summarized in Table 01 The highest values for £ have the curves 
in Figure El (b) and (c) . 
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Summary and conclusions 

In the present paper we investigated mathematical models for gene transcription 
networks with negative feedback loops. We described the time lag between tran- 
scription factor binding and mRNA production by a distributed time delay com- 
pared to previously published models which use a discrete delay. We applied our 
approach to describe experimentally observed oscillations of the transcription factor 
Hesl. 

We showed that a distributed time delay can lead to limit cycle oscillations for a 
finite range of mean delay times [f m in,T max ]. We believe this to be a more realistic 
property of gene transcription networks with negative feedback if information is 
encoded in delay-induced oscillations. More specifically, we note that any oscillation 
period should not exceed the period of fundamental cellular processes, e.g. the cell 
cycle. This is an implicit property of our model with a Hopf bifurcation and a reverse 
Hopf bifurcation leading to steady state - limit cycle - steady state transitions. 

Due to the possibility that both, sustained and damped oscillations could have a 
role in signal transmission we discussed the encoding of information in the frequency 
of the oscillations. Towards this end we proposed a quantitative measure for damped 
oscillations. 

For future experimental work we suggest a more detailed experimental study 
into the nature of oscillations in gene transcription, investigating the consequences 
of pathological mutations that lead to abnormal time lags between transcription 
factor binding and mRNA production. 
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Model 


distributed 


discrete 






r 


I: 

h p 






6 2 


6 16 


2 


8 2 


3 33 


0.2 


10 2 


2 48 


0.1 


II: 






h p 






4 2 


2 57 


0.9 


6 2 


0.6 156 


0.4 


8 2 


0.3 293 


0.2 


10 2 


0.2 470 


0.1 


6 1 


1 44 


0.4 


8 1 


0.5 91 


0.2 



Table 1: Putative Hopf bifurcation points for distributed and discrete time delays. Delay 
times are given in minutes (min). For a = 0.3min _1 and b = lmin -1 analytical results 
from eigenvalue equation Eq. © for Model I and Eq. (J2J) for Model II. We have chosen 
the degradation rate a ten-fold higher than in Ref. to ensure that the delay times are 
in the range of the time for gene transcription. 
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Model 


distributed 


discrete 




7~min 


T~max 


T 


T. 
1. 

h p 
6 2 


28 


354 


13 


8 2 


18 


558 


9 


10 2 


13 


757 


7 


10 1 


27 


163 


7 


II: 








h p 
4 2 


8 


1298 


5 


6 2 


3 


3184 


2 


8 2 


2 


5820 


1 


10 2 


1 


9208 


1 


4 1 


12 


356 


5 


6 1 


4 


991 


2 


8 1 


2 


1872 


1 


10 1 


1 


3002 


1 



Table 2: Putative Hopf bifurcation points, analytical results from eigenvalue equation 
Eq. (JSJ) for Model I and Eq. Q for Model II with a = 0.03min _1 and b = lmin" 1 . 



T 


frequency 


amplitude 


8 


0.053 


1.5 


10 


0.048 


2.0 


12 


0.043 


2.0 


14 


0.040 


1.7 



Table 3: Model I: Frequency and amplitude of limit cycle oscillations with increasing 
delay time for p = 2, h = 6 and a = 0.3min _1 . 



16 





C = w/7 




C = 


w/7 


T 


III: 


T 


I: mRNA 


I: protein 


0.3 








2 


2 


10 


2 


5 


13 


12 


50 


1 


20 


7 


8 


100 


1 


100 


3 


3 



Table 4: Frequency damping ratio £ for the curves in Figures El (Model III) and El (Model 
I)- 
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Figure 1: Model I: Intersections between the delay curves (solid lines) and the ratios curves 
(dashed lines). According to Eq. intersections indicate putative Hopf bifurcations. The 
delay curves of the distributed time delay are shown for p — 1 (in first quadrant only), 2, 
3, 4 and the ratio curves are shown for h = 2 (leftmost curve), 4, 6, 8, 10. In addition the 
unit circle as the discrete time delay curve is shown. 
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Figure 2: Model I: Numerical simulations of the exactly derived ordinary differential 
equations Eqs. (jBJ) with p = 2 and h = 6, showing a steady state - limit cycle - steady 
state transition by passing two Hopf bifurcation points with increasing mean delay time, 
(a) f = 4, (b) f = 8, (c) f = 14, (d) f = 20. Beyond the reverse Hopf bifurcation the 
damping of the oscillations increases with the delay time. For f = 50 the steady state is 
already reached at a time of 300min. 
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Figure 3: Model II: Intersections between the delay curves (solid line) and the ratio curves 
(dashed line). According to Eq. (jHJ) intersections indicate putative Hopf bifurcations. The 
delay curves are shown for p — 1 (first quadrant only), 2, 3, 4 and the ratio curves are 
shown for h = 2 (most left curve), 4, 6, 8, 10. 
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Figure 4: Model II: Numerical simulations of the exactly derived ordinary differential 
equations Eqs. (fTUj) with p = 2 and h = 4, showing a steady state - limit cycle - steady 
state transition by passing two Hopf bifurcation points with increasing mean delay time, 
(a) f = 1, (b) f = 3, (c) f = 55, (d) f = 70. 
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Figure 5: Model III: Numerical simulations of exactly derived ordinary differential equa- 
tions from Eq. (|TT|) with p = 1 and h = 4, showing no pronounced damped oscillations, 
(a) f = 0.3 (b) f = 10 (c) f = 50 (d) f = 100. Within the delay time interval [0.5,87] 
the eigenvalues are complex according to Eq. ()12j) . 
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Figure 6: Model I: Numerical simulations of the exactly derived ordinary differential 
equations Eqs. © with p = 2, h = 6 and a = 0.4, showing for intermediate delay times 
pronounced damped oscillations without passing a Hopf bifurcation, (a) f = 0, (b) f = 5, 
(c) f = 20, (d) f = 100. 
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